Boosting electronic transport in carbon nanotubes by isotopic disorder 
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The current/voltage curve of metallic carbon nanotubes (CNTs) displays at high bias a sudden 
increase of the resistivity due to the scattering of electrons with phonons having an anomalously-high 
population (hot phonons). Here, we show that it is possible to improve the electrical performances of 
metallic CNTs by ^'^C isotope enrichment. In fact, isotopic disorder creates additional channels for 
the hot-phonon deexcitation, reduces their population and, thus, the nanotube high-bias differential- 
resistance. This is an extraordinary case where disorder improves the electronic transport. 
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Introduction. — Metallic single- wall carbon nan- 
otubes (CNTs) are quasi one-dimensional wires which 
can carry the highest current density (10^ A/cm~^) d] of 
any material. This makes them the best candidates as in- 
terconnects in future high-power electronic-devices. The 
measured current/ voltage (IV) curve of metallic single- 
wall CNTs with ohmic contacts [J, 0, displays at high 
bias (voltages higher than ~0.2 Volts) a sudden increase 
of the resistivity, due to the scattering of conducting 
electrons with optical atomic vibrations (phonons) 
The largest part of this resistivity is due to the pres- 
ence of an anomalously-high optical-phonon occupation 
(hot phonons) [4|, |5|]. Indeed, the rate at which opti- 
cal phonons are excited (by the scattering with conduct- 
ing electrons) is faster than the rate at which they are 
deexcited. The large hot-phonon population determines 
an important part of the electrical resistivity, limiting 
the performances of CNTs. So far, the most direct ex- 
perimental observations of the phenomenon are given by 
Refs. [1,01, where an important increase of the optical- 
phonon population during electron transport is measured 
by Raman spectroscopy. Also, in Ref. [8] a negative dif- 
ferential resistance in tubes suspended between the two 
electrodes is interpreted as due to hot phonons. 

In this paper, we show that it is possible, in practice, 
to improve the electrical performances of metallic nan- 
otubes by ^^C isotope enrichment. Isotopic substitution 
changes the nuclear masses but not the nuclear charges 
and thus introduces elastic scattering mechanisms for 
phonons but not for electrons (i.e. it does not degrade 
the electronic transport properties) ^.9|. In fact, the in- 
troduction of isotopic disorder creates additional deex- 
citation channel for the hot phonons. This results in a 
reduction of the hot-phonon population and in an im- 
portant decrease of the nanotube high-bias differential- 
resistance. The phenomenon is quantified by means of a 
coupled Boltzmann transport equation (BTE) for both 
phononic and electronic populations, in which all the rel- 
evant scattering parameters are obtained from ab-initio 
calculations based on density functional theory (DFT). 

Interplay between electron-phonon scattering and 
thermalization. — The scattering processes involved in 



the degradation of the electronic transport are repre- 
sented in Fig.[T] In metallic CNTs, conducting electrons 
can scatter mainly with two optical phonons, correspong- 
ing in graphene to the longitudinal phonon near T 
(FEagLO) and to the near K (KA^) [13]. Their mo- 
mentum is localized in a small region of the Brillouin 
zone (BZ) [s','!^ and their energy is ^^0.2 eV. For voltages 
higher than 0.2 Volts, conducting electrons have sufficient 
energy to excite the two phonons. As a consequence, the 
electron scattering and, thus, the resistance increase. For 
tubes with a diameter of 2 nm (typically found in experi- 
ments) the rE2gL0 and KA'j^ phonons are excited with a 
characteristic electron-phonon scattering time rj^ ~ 0.5 
and ~ 0.2 ps, respectively [llj. In a defect-free 
tube, the two phonons are deexcited (or thermalized) 

with characteristic times t£, ~ 3 and rJ^ 5 ps which 

pp pp 

are due to anharmonic phonon- phonon scattering [11.] . In 
fact, anharmonicity couples the two phonons with acous- 
tic phonons which can be considered thermalized in the 
typical experimental conditions (when the tubes are ly- 
ing on a substrate) [13]. In practice, the excitation of the 
two optical phonons is faster than their thermalization 
(^ep Tpp)- Because of this, the two optical phonons 
can not be thermalized with the rest of the tube and 
their occupation is expected to increase during high-bias 
electronic transport (hot phonons). This, in turn, will 
result in a larger electron-phonon scattering and a larger 
electronic resistance. 

A quantitative description of the process can be ob- 
tained by means of a coupled BTE for both phononic 
and electronic populations, as in Ref. Q (see note [13] )■ 
The resulting IV curve reproduces very well experiments 
(compare the red curve with measurements in Fig. [2l 
upper panel). The red curve of Fig. [2] is obtained by 
letting the phonon population evolve according to BTE 
and is associated to an important increase of the opti- 
cal phonons population. To show the relevance of this, 
in Fig. [2] we also show the IV curve obtained by impos- 
ing the phonon population to be zero (blue curve). As 
already remarked in Ref. [5], in this case the high-bias 
current is much higher, that is an important part of the 
measured high-bias resistance is due to the presence of 
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hot phonons. We stress that in the calculations of Fig. [2] 
all the relevant parameters are obtained from DFT and 
no parameters are present which can be tuned in order 
to recover agreement with measurements 

The bottleneck for hot-phonon generation is the large 
phonon thermalization-time. Thus, the set up of a mech- 
anism to decrease it will diminish the hot-phonon pop- 
ulation improving the electron conduction. This goal 
could be reached by changing the isotopic composition 
of the nanotubes, since isotopic disorder provides a new 
phonon-phonon scattering mechanism. On the other 
hand, the possibility of synthesizing isotopically enriched 
nanotubes has been actually demonstrated 18, 3, 2ol |. 

In a defect-free crystal, a phonon is characterized by 
a momentum q and a branch index v. When isotopic 
disorder can be considered as a perturbation, the dy- 
namic of the system can be described by the evolution 
of the unperturbed qy states. Isotopic disorder results 
in the creation of additional scattering channels between 
the phonons, characterized by a time Tid (Fig. [1]). In 
particular, an initial hot-phonon qy can now scatter 
with numerous other phonons qV at the same energy 
hwq'u' ~ 0.2 eV (the scattering is elastic). Since such 
phonons are distributed in a large zone of the BZ, they 
are not directly coupled with conducting electrons and 
are, thus, mostly not "hot" (their occupation is thermal- 
ized and very small). In practice, isotopic disorder allows 
the scattering of a hot phonon with other cold phonons 
adding a deexcitation channel and, thus, increasing the 
thermalization rate from t^^ to t^^ + . 

Computation of the isotopic- disorder scattering- 
times. — Let us consider an isotopically enriched tube, 
with disorder parameter x. m (m-f Am) is the mass if the 
i2(j (^^C) isotope. The atomic mass nii is m-\- Am with 
probability x, being in average m = m -t- xAm. In per- 
turbation theory with respect to Am/m, the lifetime 
due to isotopic disorder of a phonon with pulsation 
Wq^i/ is given by the Fermi golden-rule 



21 



qi/N_l 



52 ( 



(1) 



where the sum is performed on Nq q-points, s is an 
atomic index, eqy{s) is the phonon polarization nor- 
malized as |eqy(s)p = 1 in the unit-cell and 5,, is a 
Lorentzian with full-width at half maximum 77, which is 
determined self-consistently from 77 = [t^^ iini)]~^ ■ g2{x) 
describes the mass disorder and is 
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and Gqi/ are computed from DFT |12l. |23| 



To verify the validity of Eq. [1] t'^^ are also determined 
exactly, following Ref. 25]. Let us consider a super-cell 



containing N atoms and let us call D = M~^/^KM~^/^, 
where M is the diagonal N x N matrix describing the 
random distribution of the masses (M^- = rriiSij) and 
K are the interatomic harmonic force constants. The 
phonon spectral-function is given by: 
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where V^y are the atom displacements of the unper- 
turbed qz/ phonon (normalized to 1 in the super-cell) 
and ^ is a small positive number. Using the DFT force- 
constants [12, 23| and computing A(qi/, w) with a Lane- 
zos method [2 

full- width at half maximum of A[(\v, lo 



T^^ are obtained as the inverse of the 
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Tid are computed for several tubes with different diam- 
eters d with both Eqs.[T]and [3l Fig. [3] shows the results 
for the two phonon branches relevant in electronic trans- 
port for a rf = 2 nm tube. The results from the two 
methods are very similar. Thus, the perturbative ap- 
proach (Eq. [1]) is a good approximation to the exact one 
(Eq. [3]). Moreover, the results obtained for a d = 2 nm 
nanotube are very similar to those for two-dimensional 
graphene. Thus, the results for tubes with d > 2 nm are 
not strongly dependent on the diameter. 

It is now important to remark that the Tid for the two 
hot-phonons (r^^ and t^ in Fig. [T|) are smaller than the 
corresponding Tpp. Thus, the new thermalization chan- 
nels (due to isotopic disorder) are more efficient than 
those present in a defect-free tube (due to anharmonic- 
ity) and are expected to influence the IV curve. 

Computation of the IV curve. — To quantify the im- 
pact of isotopic disorder (ID) on the high-bias electrical 
resistivity, we use the BTE as in Ref. [5|]. The ID scat- 
tering rate of an initial qy phonon can be decomposed as 
a sum of scattering rates into the other qV phonons, by 
rewriting Eq. [T] as (rj)"! = J^q',.' V^-,qV The ther- 
malization process due to ID is then taken into account 
by adding to the collision term for the phonons (Eq. (4) of 

Ref. i) [dtn'^mso = Eqv' r^Lq'.'i^^'i^') " ^"^(q)), 
where n''(q) is the occupation of the qv phonon. 

First, Fig. [4] reports as example the ID transition rates 
for a particular hot-phonon qv. The majority of the qV 
states are not "hot" (we checked that more than the 90% 
of the scattering rate is due to phonons which are not 
"hot"). Second, for each of the qV phonons, Tqi^^qi^i 
is much larger than their anharmonic scattering rate (as 
obtained in Ref. (lj|). A consequence of these two facts is 
that the population of the phonons which are not directly 
coupled to conducting electrons, n'^ (q'), is negligible and 
one can rewrite [dtn''{q)]^^^ ~ -(Eq'i.' V^-^pi^O'^''(q)- 
That is, one can directly use Eqs.1-4 of Ref. [5|, replace 
the thermalization time Tj^^ by {T^p)~^ + {t^)~^ cor- 
responding to the two hot-phonons and consider the q 
dependence of Tid [IBl ■ 

In Fig. [2] we report the IV curve for the isotopically en- 
riched (a;=0.5) carbon nanotube. Indeed, the conduction 
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properties are ameliorated with respect to the isotopi- 
cally pure case {x = 0). To be quantitative, we notice 
that at high-bias the differential resistance dV/dl is al- 
most independent from the voltage and is linear with 
respect to the tube length (lower panel of Fig. [2]). This 
allows to define a resistivity p and an effective scattering 
length Isc, by 
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where ft./(2e^) — 12.9 kil is the quantum of resistance. 
From Fig. [21 the introduction of isotopic disorder en- 
hances the scattering length ( and thus reduces the high- 
bias differential resistance) by a factor three. 

Conclusions. — The introduction of isotopic disorder 
in metallic carbon nanotubes creates an effective channel 
for the thermalization of the hot phonons without intro- 
ducing an additional scattering process for the conduct- 
ing electrons. This is expected to reduce the actual elec- 
tronic differential-resistance. In particular, according to 
a Boltzmann transport treatment (in which all the rele- 
vant parameters are obtained from ab-initio and no semi- 
empirical parameters are present which could be tuned 
to reproduce the experimental results), isotopic disorder 
reduces the high-bias electronic differential-resistance by 
a factor of three in tubes with a diameter of 2 nm. This 
amelioration of the transport properties is expected to 
have important technological consequences in view of 
the use of metallic carbon nanotubes as interconnects 
in tomorrow electronics. Further improvement could be 
reached by adding other sources of phonon scattering 
wich arc transparent to electrons as, e.g., by cycload- 
dition functionalization [2^. 

After submission of the present work, another paper 
appeared proposing isotopic enrichment to reduce hot- 
phonon population in GaN transistors |28j]. We thank 
G. Stolz, N. Bonini and N. Marzari for discussions and 
aknowledge IDRIS project 081202. 
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FIG. 1: Schematic view of the hot phonon generation (see the 
text). In a defect-free tube, the hot phonons are generated 
and thermalized with characterisitc times Tep and Tpp, respec- 
tively. Isotopic disorder introduces additional thermalization 
channels (shaded area) with characteristic times nd- 
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FIG. 2: (Color online) Upper panel: IV curve of a (14,14) 
nanotube with length L and diameter d. Black lines are mea- 
surements from Refs. ^ |3| . Red and violet lines are calcu- 
lations for an isotopically pure {x = 0) and enriched tube 
{x = 0.5). The blue line is obtained by imposing the phonon 
populations to be zero. Lower panel: high-bias differential 
resistance for (14,14) tubes with different lengths (same color 
code as in the upper panel). Lines are linear fit to the points 
and Isc is the resulting scattering length (Eq. [4]). 



do not depend on the nuclear masses. Thus, isotopic sub- 
stitution scatters electrons only by mechanisms (which in 
general are negligible) beyond the Born-Oppenheimer ap- 
proximation and/or beyond the harmonic approximation 
for interatomic force-constants and for electron-phonon 
couplings. These mechanisms have been studied in liter- 
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FIG. 3: Phonon lifetimes due to isotopic disorder (a; = 0.5) in 
a (14,14) CNT and in graphene computed with the perturba- 
tive (pert.) and exact approach of Eqs.[T]and |3l respectively. 
The phonons in the two panels correspond to the highest op- 
tical branches of graphene including the rE2gLO and KA'^ 
modes. The grey areas correspond to phonons with high oc- 
cupation (hot phonons) at 1 Volt. 




FIG. 4: (Color online) Transition rates T^J'_qi^i from an ini- 
tial qy phonon (white dot) into qV phonons (see the text) 
due to isotopic disorder {x = 0.5). Horizontal lines are the 
wavevector in graphene Brillouin zone (BZ) compatible with 
the (14,14) CNT. Filled curves are proportional to t~J'_^^,^i 
The red zones correspond to phonons with high occupation 
(hot phonons) at 1 Volt. 
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dicted electron scattering-lengths of the order of the mil- 
limeter. These scattering lengths are undetectable with 
respect to those due to other sources of electron scatter- 
ing found in bulk metals and in carbon nanotubes, where 
the relevant lenghts range from 10 nm (electron-phonon 
scattering) to 1600 nm (elastic scattering) 
[10] In metallic CNTs, the electronic bands cross the Fermi 
energy at two equivalent points of the reciprocal space 
corresponding to the symmetry points K and K' of the 
graphene BZ [j . Because of momentum conservation, the 
conducting electrons of the tube can scatter only with 
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